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We discuss the high energy neutrino emission from gamma-ray bursts resulting from the earliest 
generation ('population III') stars forming in the Universe, whose core collapses into a black hole. 
These gamma-ray bursts are expected to produce a highly relativistic, magnetically dominated jet, 
where protons can be accelerated to ultra-high energies. These interact with the photons produced 
by the jet, leading to ultra-high energy photo-meson neutrinos as well as secondary leptons and 
photons. The photon luminosity and the shock properties, and thus the neutrino spectrum, depend 
on the mass of the black holes as well as on the density of the surrounding external gas. We 
calculate the individual source neutrino spectral fluxes and the expected diffuse neutrino flux for 
various source parameters and evolution scenarios. Both the individual and diffuse signals appear 
detectable in the 1-300 PeV range with current and planned neutrino detectors such as IceCube 
and ARIANNA, provided the black hole mass is in excess of 30-100 solar masses. This provides a 
possible test for the debated mass of the progenitor stellar objects, as well as a probe for the early 
cosmological environment and the formation rate of the earliest structures. 

PACS numbers: 



I. INTRODUCTION 

The first generation of stars in the Universe (known 
as population III stars) are the earliest objects to form 
from the collapse of pristine gas. There is so far no di- 
rect observational evidence for their existence or prop- 
erties, but simulations [1-4J suggest that they are likely 
to form with a top-heavy initial mass distribution, with 
a much heavier average mass than current stars, in the 
range 10^ — IQ^Mq, although models with smaller masses 
are also possible, e.g. [S]. Such hypermassive stars have 
very short life times, and except for a limited intermedi- 
ate mass range, they undergo a core collapse leading to 
a black hole whose mass is some fraction of the initial 
stellar mass, depending on how much mass loss occurred 
during the brief stellar evolution. For sufficiently fast ro- 
tating cores, accretion of remnant gas onto the black hole 
can lead to the formation of a powerful jet, resulting in 
a gamma-ray bursts (GRB), e.g. [SHE]. 

In this paper, we discuss a specific scenario where the 
GRB jet is initially dominated by magnetic energy (an 
MHD jet) [TU]. A very large energy output can be 
realized in the jet, leading to shocks which can accel- 
erate protons to highly relativistic energies, leading to 
significant neutrino emission via the photomeson process 
in the presence of the accompanying large photon lu- 
minosity. Although difficult to observe, these neutrinos 
propagate almost absorption-free across cosmological dis- 
tances, and can provide valuable information about cos- 
mic conditions during the reionization epoch, when the 
first structures formed. 

In this paper we discuss the very high energy neutrino 
flux expected from Pop. Ill GRBs and the potential for 
its detectability with current or future large neutrino de- 



tectors such as IceCube [II], ANITA [H] and ARIANNA 
[T^ . In section |ll] we discuss the method of calculation. 
The main features for our Pop III GRB model and its as- 
trophysical setting are described in subsection Ila. The 
photon spectrum serving as the target for the photome- 
son process is discussed in subsection lib. The method 
used to calculate the photomeson process and other com- 
peting channels for proton and secondary particle energy 
losses is detailed in subsection lie. In section III we dis- 
cuss the possible Pop III GRB source rates, and the po- 
tential for the detection of the corresponding neutrino 
fluxes. A summary and discussion of the implications is 
given in section IV. 



II. MODEL AND CALCULATION METHOD 

A. Astrophysical Input 

In the very massive 10^ — IQ^Mq Population III star 
scenario, aside from the approximate 140 — 26OM0 range 
leading to pair instability supernovae which leave no rem- 
nant [7| , the core collapse results in a massive black hole 
(BH) of mass Mh encompassing a substantial portion of 
the original mass. A fraction of these are expected to be 
fast rotating ^141 [l5], one of the requirements for GRB 
progenitors. For a star of initial radius i?*, infall of the 
remnant stellar gas onto the black hole leads to an accre- 
tion disk of typical outer radius Rd ^ -R* /4 which for a 
typical disk magnetization parameter a = 10~^a_i gets 
swallowed on a timescale |5j 



2 



tn 



7 ( m 

3a ^ GMh 



1/2 



1.4 X 



10^aZlRl%M- 



1/2, 
2.5 " 



(1) 



where Af/j 2.5 = 10 is the largest mass of the cen- 

tral black hole considered here, and henceforth we use 
the convention Ax = . Extraction of the rota- 

tional energy from high mass fast-rotating black holes 
is expected predominantly through MHD effects via the 
Blandford-Znajek mechanism jl6j . The corresponding 
BZ luminosity of the resulting jet is estimated as 



L 



BZ 



128 



(2) 



Here is the dimensionless spin parameter of a Kerr 
BH , Rh ^ GMhjf? and Bh is the disk magnetic field 
strength threading the black hole horizon, which should 
scale with the disk gas pressure P in the advection- 
dominated disk (ADAF) regime [l2], = where 
/3 = 10 Thus, Bh ~ [A^/liMc/iaPRlf^ ~ 6.6 x 
1013 {a^iPi)~^M-}^^My^l^t~]/'^ G, where is the ac- 
cretion disk mass. Assuming a simple scaling Md = SMh, 
where S ^ 1, the total energy output of the jet is 



E, ~ LBztd ^ 2.2 X 105^(af,(5/a_i/?i)M,,,2.5 erg, (3) 



and the isotropic equivalent energy is 

E,so = E,{1 - cos0j)~i ~ E,i2/9' 



(4) 



where (1 — cos6'j) ^ ^|/2 is the beaming factor for a jet 
of coUimation half- angle 6j ^ 1. 

The ratio of the photon scattering mean free path and 
the jet radial dimension in the comoving frame (the scat- 
tering "optically depth') is very large near the black hole 
and drops outward, until it becomes unity at a photo- 
sphere located typically beyond the original stellar ra- 
dius. As the jet expands beyond this photospheric ra- 
dius, if the jet is strongly magnetically dominated, in- 
ternal shocks are unlikely to occur [9l [10], but an an 
external forward shock (FS) forms as the jet sweeps up 
the external gaseous matter in its surrounding. For an 
approximately uniform external density of n/cm^, at the 
time td when the accretion process stops feeding the jet, 
the jet head has reached a distance Vd from the central 
explosion. 
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and the bulk Lorcntz factor of the jet head is 
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e.g. [H], where from now on we will write E = E^go, 
and nip is the proton mass. In the standard GRB shock 



description, the random magnetic field energy density in 
the comoving frame of the shocked external gas is ampli- 
fied to some fraction of the internal energy, 

B ~ {^2^tBnmpC^f'^Td ~ 2,.s4^_^El!,%t-/J\l'^G. 

(7) 

In the shock the electrons are Fermi-accelerated into 
a power law distribution, resulting in synchrotron emis- 
sion in the above field, and these synchrotron photons 
are further subjected to scattering by the same electrons 
leading to a synchrotron-self-Compton (SSC) radiation 
field [lU]. 

If by the time td the original jet magnetization parame- 
ter has decreased, which can also be promoted by baryon 
entrainment from its surroundings, besides the forward 
shock (FS) also a reverse shock (RS) will develop, which 
moves into the ejecta |18j . and a hydrodynamical approx- 
imation can be used for the description. Let the Lorentz 
factor of the unshocked jet in the source frame be F^-, 
the jet head Lorentz factor be F^ (also measured in the 
source frame) and the Lorentz factor of the unshocked 
ejecta measured in the jet head frame be F*. Since in 
our case Fj , F^ ^ 1 , we have 
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The comoving number density in the unshocked ejecta is 
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The Lorentz factor of the jet as inferred from observed 
GRB afterglows is 10^ ^ IQ'^, and in these calculations 
we adopt Tj 500 as a nominal value. The majority 
of the reverse shock emission is produced when the re- 
verse shock finishes crossing the ejecta and injection of 
fresh electrons ceases [19,. This shock crossing time is 
essentially td, at which time the ejecta thickness is ap- 
proximately ctd and the shock radius is approximately 

The details of the photon spectrum vary in time and 
depend on the energy of the GRB jet, the shock physics 
parameters and the details of the environment. The 
fraction of the shock energy that goes into relativistic 
electrons Ce (the electron equipartition parameter) is as- 
sumed to be the same for the forward and reverse shocks; 
a similar assumption is made for the forward and reverse 
magnetic field and the accelerated proton energy equipar- 
tition parameters es and Sp (relevant for neutrino emis- 
sion). 

The jet total energy Ej is reasonably well defined as 
a function of the disk mass and the BH spin parameter. 
However, the corresponding isotropic-equivalent energy 
Eiso depends on the uncertain jet opening angle. Fi- 
nally, the radiation produced depends, via the shock ra- 
dius, on the external medium density n, which is largely 
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unknown, depending on details of the star formation pro- 
cess. For simplicity we treat this density as a constant 
free parameter, our choices being guided by existing nu- 
merical simulations of early star formation and the re- 
sults of analyses and modeling of lower redshift GRBs. 
The modeling of observed GRB afterglows suggests that 
for the z < 8, i.e. later generation, so-called Pop. I/II 
GRBs, the densities of the medium in their environment 
typically range over 0.1 < n < 100 cm~^ ISUj. The envi- 
ronments of the first stars prior to their collapse has so 
far only been inferred from model numerical simulations, 
which differ significantly among each other. For example, 
the typical early galactic gas environment could evolve 
as n cx (1 + z)' |21j . or it might be approximately inde- 
pendent of redshift, n ~ 0.1 cm~'^, as a result of stellar 
radiation feedback [521 US] ■ 

The small number of analyses for what are currently 
the most distant GRBs imply ambient densities for these 
high redshift GRBs which could be n ~ 10^ — lO'^cm"'^ 
for GRB 050904 at z ~ 6.3 [H], and n ~ Icm^^ for GRB 
090423 at z ~ 8.2 [25]. 

As nominal cases, we will discuss mainly the models 
in Table 1, the most energetic examples having — 
3OOM0, {al6/a-i(3i) ~ 1, Ej = lO^^ ^ erg, es = 0.01, 
Ep = 0.1, jet opening angles 6j = 10~^, 10^^ and an 
external density n — 1, 10^, 10' cm^'^, the corresponding 
isotropic equivalent energies being given by eq. Q. We 
also consider smaller black hole masses Mh — IOOA/0 
and 30Mq with correspondingly lower isotropic energies. 
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Table 1 



Table 1 caption: Pop. Ill GRB model parameters 



These parameter values are reasonable, given the various 
uncertainties, but are by no means unique. Different en- 
ergy equipartition parameters ex are possible, and also 
lower Pop. Ill stellar masses [5] , which would lead to even 
lower Mh and Ej than those in the table. Such lower val- 
ues, however, will produce dimmer neutrino fluxes (see 
below) which are not favorable for detection with cur- 
rent or future neutrino detectors. Higher Mh (e.g. [26] ) 
and Ej than those in Table 1 are also possible but more 
speculative. 



B. Photon spectra 

The GRB photons provide the most abundant targets 
for the photomeson process (nuclear collisions, despite 
larger cross sections, are much rarer). For times t < 
td, a common feature of all jet models is a photospheric 
spectral component, arising at a radius Vph much smaller 
than r^, where the photon scattering timescale equals 
the expansion timescale. This component has a quasi- 
blackbody spectrum[101, and its fluence is taken to be a 
fraction =0.1 of the jet total energy. 

In addition to this quasi-thermal component, the 
shocks contribute various non-thermal photon spectral 
components, which typically dominate at energies both 
above and below the photospheric component. Electrons 
are accelerated in the shocks to a power law spectrum 
(in/d'ye for 7e > 7m, where 
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is the minimum injected electron Lorentz factor in the 
shock comoving frame (where i = *, c? is for reverse and 
forward shock, respectively), and f{p) ~ (p — 2)/(p — 1) 
since observations suggest an index 2 <p < 2.5. Acceler- 
ated electrons produce a (comoving frame) synchrotron 
spectrum peaking at 
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and the synchrotron photons are inverse Compton (IC) 
scattered by the same electrons to produce a synchro-self- 
Compton (SSC) spectrum peaking at 
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The different parameters of the cases A-D in Table 1 
lead to non-thermal photon spectra which can be quite 
different. In case A, the non-thermal photon spectrum 
consists of the sum of the synchrotron emission and the 
SSC component. In the cases B-D the energy density of 
the photons is significantly higher than in case A. The 
peak of the original SSC spectrum exceeds the threshold 
energy for two-photon pair formation, E^ ^ E^^. Here 
Ejj is the 77 self-absorption energy above which high en- 
ergy photons produce cascades of electron-positron pairs. 
The secondary pairs in turn also produce synchrotron 
emission and IC-scatter photons. Part of the SSC spec- 
trum, however, can be suppressed by the Klein-Nishina 
effect. 

In cases B and C, the effects from one generation of 
cascade pairs are calculated. The results suggest that 
the second generation of pairs affects the spectrum less 
significantly than the first and is neglected for simplic- 
ity. In case D the photon energy density is so high in 
the higher mass sub-cases that the copious e* created by 
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e can modify the original photon spectrum by 



multiple Compton scatterings. The compactness param- 
eter (defined as the comoving optical depth to 77 effects) 
can be estimated as 



(13) 



for the case I?3oo> and about I' ^ 34 for I?ioo, so the 
spectrum is partly thermalized (it would be completely 
thermalized if I' ^ 10^, e.g. [57] ). In addition, hadronic 
cascade photons from neutral pion decay could also af- 
fect the spectrum, an exact evaluation of such cascades 
requiring a numerical calculation |27| . In general, how- 
ever, the effects of hadronic cascades is sub-dominant 
relative to electromagnetic cascades from 77 effects [28]. 
To simplify things, in this paper we calculate the photon 
spectrum of case D under two limits: -D(l), where the 
non-thermal spectrum is approximated as a summation 
of the synchrotron and SSC spectrum, and D{2), where it 
is assumed that cascades lead to a completely thermal- 
ized (black body) spectrum. Due to the uncertainties 
in case D, only the forward shock emission is considered. 
(As it turns out, the neutrino spectra calculated for cases 
D{1) and D{2) give similar predictions for the flux in the 
ARIANNA energy band of observational interest, since 
in this energy range the pionization efhciency approaches 
unity; see SIII). 



The details of these various photon spectra are cal- 
culated using the methods described in the appendix of 
|10j . The resulting photon spectra for all cases in Table 
1 were calculated. In Figures 1-3 we show only some 
of the representative cases A, B and D, expressed in 
the observer frame, for the mass sub-case Mh = 300Mo. 
The smaller mass cases are roughly similar, downscaled 
versions of these. 



C. Proton Acceleration, Cooling 
and Photomeson Production 



density in the comoving frame. The total inverse Comp- 
ton time scale is given by j3lj . 
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where F{E, 7) 

f2{zb)], Za = 

{E/'mpc'^){'y+y/ — 1)""^, and fi{z), f2{z) are functions 
given in f3T|. The rclativistic protons interact mainly 
with GRB photons, which are the most abundant tar- 
gets in the acceleration region, including photospheric 
photons as well as forward and reverse shock photons. 
When the feeding stops at t ^ the photospheric and 
reverse shock emission decay very rapidly, (unless there is 
continued central engine outflow), and the forward shock 
decays as power law afterglow. Thus, the bulk of the pho- 
ton production and proton acceleration peaks around the 
time td, and here we concentrate primarily on the emis- 
sion over this timescale. 

There are various channels for p7 interaction, the 
two most important ones in our case being the photo- 
meson and the photo-pair channels. In the photo-meson 
channel, neutrinos are mainly produced by the decay of 
charged pions, Kaons, etc. For the two nominal cases 
considered in this paper, the ultra-high energy protons 
have a low probability to interact with high energy pho- 
tons, and Kaon production is not significant [0] . 

In the photo-meson mechanism approximately equal 
numbers of charged and neutral pions are produced, 
p + 7 — > j3 + 7r° and p + 7 — )• n + 7r+. The neutral 
pions decay into very high energy photons, while the 
charged pions lead to neutrinos tt"*" — >■ /i+ -I- with 
/i+ e'^ + Vf, + i^fj.- In this paper, we use a two-step- 
function approximation for the photo-pion production 
cross section and for the averaged inelasticity to approx- 
imate single-pion production and multi-pion production. 
In this approximation the cross section is 



Protons will be Fermi accelerated at the shocks to 



form a spectrum N cc E e.g. 



The acceleration 



timescale for protons can be estimated as 
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where ~ 1 — 10 is a shock-structure dependent coef- 
ficient [30]. 

The adiabatic cooling time scale is tadb ~ ^d/Tc . The 
accelerated proton is also subject to synchrotron cooling 
with a time scale 



f340Aib for ethr = 390 < < 980 
'"^^^'^^-j 120Mb for ..> 980 ' ^^^^ 

where is the invariant energy, and the inelasticity is 



^ , , J 0.2 for 390 < < 980 ^, 
0.6 fore. > 980 ^^'^ 



tsyn SfTlpC /AcTTm^EpUB 
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where ax = 0.664/xbarn is the Thomson cross section 
for electrons and Ub — B'^/Stt is the magnetic energy 



(although there may be a small pile up at GZK energies (where 
charged pions suffer significant radiative cooling; while Kaons 
suffer less), which is beyond the scope of the present calcula- 
tion and would require a Monte Carlo simulation toolkit such as 
GEANT4 [62]). 
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Figure[l]: Photon E^{dN/dE) spectrum in the ob- 
server frame from a Pop. Ill GRB at z = 20 , case 
^300 (Figs. 1 through 6 are for Mh = SOOM©). SoHd 
Unes- left: synchrotron component, right: SSC compo- 
nent, both from the forward shock. Dotted: photospheric 
emission. Solid envelope: total of forward shock plus 
photospheric emission. Dashed lines- left: synchrotron 
component, right: SSC component, both from the re- 
verse shock. Dashed envelope: total of the reverse shock 
emission. 
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Figure[2]: Photon E'^{dN/dE) spectrum in the ob- 
server frame from a Pop. Ill GRB at z = 20 , case 
i?30o- Solid lines- left, syn(FS): synchrotron compo- 
nent, right, SSC(FS): SSC component, both from for- 
ward shock. Dotted: photospheric emission. Other solid 
lines labeled "p" : emission from electron-positron pair 
cascades. Overall solid envelope: total forward shock 
emission, included photospheric emission and cascade. 
Dashed lines- left: synchrotron component, right: SSC 
component, both from reverse shock. Dashed envelope: 
total reverse shock emission. Thick vertical dashed line: 
photon-photon absorption energy. SSC spectrum (FS) 
above E^i^i (Vertical dot-dashed line) is suppressed by 
the Klein-Nishina effect. 
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Figure [3]: Photon E'^{dN/dE) spectrum in the ob- 
server frame of a Pop. Ill GRB at z = 20 , case I?3oo- 
Two extreme limits are shown: Case D(2) assumes com- 
plete thermalization of the original spectrum due to copi- 
ous pair formation and Compton scatterings; Case D(l) 
solid lines show the original synchrotron and SSC from 
the forward shock, dashed is the photospheric emission, 
and the solid envelope is the sum of these, neglecting 
spectral changes due to pair formation. 



The detailed prescription and the ratio of these prod- 
ucts are described in [52] ■ The photohadronic cooling 
time scale is 



27p 
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(19) 

where dn/de is the photon number density per energy 
decade in the shock comoving frame. 

The other important pj interaction mechanism under 
the present conditions is photo-pair production, p + 7 — > 
p+e'^+e~ . Although in general the energy loss rate from 
this channel is small compared to photo-pion production, 
due to the fact that the average inelasticity per collision 
is small ( ^ 10~^ ), in some cases and at some energies 
multiple interactions can take place and the energy loss 
rate can become appreciable compared to the photo-pion 
rate. Here we treat this effect using an approximate an- 
alytical method described in [33j, taking the photo-pair 
cooling time scale to be 
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where^ = {TnieafcaT I^V^T^rnp) is a constant, a/ ~ 
1/137 is the fine structure constant and /c^e = 1 is a 
numerical constant selected to approximate the relevant 
cross-section. Charged pions and muons are also subject 
to energy loss due to radiative cooling. The dominant 
channel is the synchrotron cooling. Due to this effect the 
logarithmic slope of the neutrino spectrum steepens by 
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2 units [ini [30] above the muon critical energy E'^j, at 
which 

tsy,niE^b) = tdec,ti{Ef^b), (21) 

where t^y^^ is the muon synchrotron cooHng time scale 
and tdec,ij. is the muon decay time scale in the comoving 
frame. The neutrino spectrum is further suppressed 
above the pion critical energy E^^f, where £^7r6 is given 
by tsy,-K{ET,i,) = tdecay,AE-Kb)- Rclativistic neutrons 
produced through the channel p + 7 — > n + 7r+ have a 
much longer decay timescale [34] in the comoving frame 
and their effect is neglected in this calculation. 

The re-acceleration timescale for pions is longer 
than the pion decay timescale for all cases of interest, 
while for muons t^^ acc < tij.,decay when the magnetic field 
is above ^ 5 G . The latter can be realized in the early 
phase of the afterglow in our B and D case. However, 
for simplicity, we shall assume here that re-acceleration 
is inefficient for all leptons when calculating the photon 
and neutrino spectrum. A more detailed spectrum in- 
cluding the re-acceleration can be evaluated numerically, 
as in [19]. 

The total energy loss rate of the protons is given 
by = '^t~^{i = all channels) and the photo-pion 
cooling efficiency is defined through fnilp) = tp^ /ip^- 
This gives the average fraction of energy lost to pions 
from the injected protons at energy 7^. On average each 
charged pion carries a fraction ^ 0.2 of the energy of 
its parent proton and each neutrino carries a fraction of 
^ 0.05 (either from pion decay or muon decay). Thus 
we have 

J,{E,) - {l/^)h{Ep)U.^b{Ep)Jp{Ej,) (22) 

where E^ ^ Q.QbEp and the flux Jx is dcflned by 
Jx = E\dN{Ex)/dExdt . The function UM^p) - 
Mtn[l, iiEjE^b)-^)]{[il/2)Mtn[l, (E^/E^b)-^] + 
(l/2)Afm[l, ((i?^/i?^fc)^^)]} approximates the effects 
of pion and muon cooling discussed above. The use of 
the factor 1/2 assumes that half the neutrinos come 
from charged pion decay and the other half from muon 
decay. Neutrinos with different flavors approximately 
contribute equally when they oscillate [34] over cosmo- 
logical distances although muons are more cooled than 
pions before they decay and may induce a different 
flavor ratio 1351 



In flgures 4-6 we show the comoving frame cooling 
timescales for the various proton interaction channels, 
corresponding to photon spectra of Figs. 1-3, for the 
forward shock regions in the cases A^qq, S300 and 

i^(2)300. 
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Figure [4]: Proton cooling, case A300 inverse timescale in 
the forward shock region, plotted in the jet comoving 
frame. The thickness of the acceleration line shows the 
uncertainty in acceleration efficiency. 
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Figure[5]: Proton cooling, case i?3oo inverse timescale 
in the forward shock region plotted in the jet comoving 
frame. 
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Figure[6]: Proton cooling, case 15(2)300 inverse timescale 
in the forward shock region, in the jet comoving frame. 
The cutoff in the dashed line is due to the photomeson 
threshold and the cutoff in the dotted line are due to the 
photo-pair threshold, where the photon energy in the 
proton frame are ~ ^ 140MeV and ^ ■m^± ^ IMeV 
respectively. 
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III. NEUTRINO FLUXES 



Individual source neutrino fluence 



Pop. Ill stars are expected in the re-ionization epoch, 
at redshifts z'>7 [3S1I37], and the majority could arise 
at redshifts z ~ 20, with the very first objects at red- 
shifts possibly as high as z < 70 [SS]. Here we calcu- 
late the neutrino spectrum for an individual Pop. Ill 
GRB using the model and approximations described in 
sections I and II, for the various cases listed in Table 
1 of section II. The total neutrino fluence spectra (the 
time integrated energy spectral flux) are evaluated at 
the deceleration time td, when the GRB emissivity is 
the largest. As an example, the source fluences for the 
highest mass cases Mh — 3OOAf0 are plotted in Figure 
7 in the observer frame, for the four nominal cases A300 
through _D(1, 2)300, assuming a redshift z = 20. At low 
energies the spectrum is dominated by the A-resonance 
where the proton energy reaches the threshold energy. At 
higher energies multi-pion processes become significant. 
The total emission is dominated by the interaction with 
the external shock photons, especially synchrotron pho- 
tons (due to the original electrons and, depending on the 
case, secondary pairs), for which the photo-pion efficiency 
approaches order of unity ~ o(l) . The IC and SSC pho- 
tons contribute less significantly due to their low number 
densities. This is the reason why the shape of the total 
neutrino spectrum is similar between the cases A — C . 
However in case D{2), where photons are assumed to be 
complete thermalized, there is a dearth of high energy 
photons above the Wien tail so only very high energy 
protons meet the photo-pion threshold, and the neutrino 
spectrum is different from those in cases A — C. 



As seen in Figure 7, the neutrino emission from such 
Population HI MHD-dominated GRB models is charac- 
terized by a very high peak flux energy, as well as a 
high fluence. The latter is due to the high black hole 
mass, which implies a high intrinsic luminosity as well 
as a long duration of the external shock peak emissivity 
phase (which is further lengthened in the observer frame 
by the high redshift). The magnetic field strength at 
the external shock is weak compared to that in internal 
shocks, so the cooling effect for charged pions becomes 
much less significant than when internal shocks may be 
important, such as in lower redshift GRBs. An inter- 
nal shock component is not included here, since internal 
shocks are unlikely in strongly MHD dominated GRBs, 
|2] (for internal shock neutrino emission in hydrody- 
namical models see, e.g. [Tnil3SH12)- 




Log,„|£;.l (GeV) 

Figure [7]: Neutrino fluence (time-integrated energy flux) 
from one Pop. Ill GRB of Mh = SOOM© located at 
2; = 20 . Model cases ^300 through C300 are the labeled 
solid lines, model D is shown as dashed lines for the 
two extreme cases D{1) and D{2) corresponding to the 
photon spectra of Fig. 3. The neutrino emission from 
the forward shocks is shown by thick lines and that from 
the reverse shocks as thin lines. The thick dashed line 
is the neutrino emission in case A (forward shock) if 
photospheric photon emission were absent, to show the 
effect of the latter. Cases D(l) and D(2) result in very 
similar fluence levels in the multi-PeV range, because 
the efficiency for photo-meson pion production is nearly 
unity at these energies. 





E 


B 


D{2) 


D{1) 


10 1 20 1 70 


10 1 20 1 70 


10 1 20 1 70 


300 


TeV - PeV 
PeV - EeV 


3.0 2.1 1.4 

6.1 3.2 1.0 



2.7 1.5 0.65 


10.3 5.6 2.1 
3.8 1.6 0.37 


100 


TeV - PeV 
PeV - EeV 


0.71 0.51 0.33 
1.6 0.83 0.28 


0.64 - 


4.5 2.7 1.1 
2.1 0.90 0.18 


30 


TeV - PeV 
PeV - EeV 


0.16 0.11 - 
0.41 0.22 - 




1.8 1.0 0.40 
0.72 0.30 0.06 



Table 2 Caption: Number of muon events from an 
individual burst, in models B and D for ~ 
300, 100, 30Mq from redshifts z = 10, 20, 70. 

We can obtain an initial quick estimate of the number 
of muon events expected in a km^ detector from one burst 
such as, e.g., cases B or D{1), D{2) in Fig. 7. Taking 
the fluence to be E^{dN/dE) - 1 GeV cm^^ in aU fla- 
vors in a band AE - 10^-^ GeV around E - lO^-^ GeV, 
there are AA^i, ~ 3 x 10^^ cm~^ neutrinos received over 
a time td^obs ^ 3 x 10^ s, or A^j, ~ 3 x 10^ km~^ per burst. 
Complete mixing occurs over cosniological distances, so 
1/3 of those are muon neutrinos, and for an approxi- 
mate conversion probability P^-j.^ ^ 1.3 x 10~^{E^/TeV) 
(valid in the range TeV Ei, < PeV) this translates 
into AN^ ^ 4 km^'^ muon events per burst from a 
redshift z = 20. A more accurate calculation takes 
into account that the conversion probability changes to 
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P^^f, ~ 10-2(i;^,/EeV)0 '*^ in the range PeV < < 
EeV, where EeV= 10^^ eV. Table 2 shows the result of 
such a calculation performed numerically using the de- 
tailed spectra of Fig. 7, which gives the number of muon 
events separately in the TeV-PeV and PeV-EeV energy 
ranges. 

The above are only the muon events, in addition to 
which there would be tau events, which can increase 
the signal (up to at most a factor 2 at the highest en- 
ergies). These Pop. Ill GRB neutrino signals need to 
be considered against at least three different sources of 
background, (a) One is the diffuse atmospheric neutrino 
background, which has a steep spectrum, and at energies 
~ lO^-S GeV has an upper limit of gS] E'^^e ~ 10"^ 
GeV cm~^ s~^ sr"-*^. Taking an angular resolution circle 
of ~ 0.7 degree, this gives over the duration of the burst 
AiV,,atm 6 X 10-5 km-2 or A7V^,atm - 2.55 X 10-^ 
km~^ muons per burst, a negligible background, (b) 
Another background is the GZK cosmogenic neutrino 
background, due to the photo-meson interactions of the 
observed ultra-high energy cosmic rays with the cosmic 
microwave background photons [1^49j . While the ex- 
act value depends on the assumed evolution with red- 
shift of the cosmic ray sources, it is generally important 
only at energies higher than those where the signals pre- 
dicted here could be important, (c) A third background 
is that which may be expected from lower redshift (Pop. 
I/II) GRBs, the typical model for which [551 [30] is cur- 
rently close to being constrained by IceCube measure- 
ments [32 ED. 

From Table 2 one sees that for masses Mh ^ 3OOAf0 
at z = 20 the signal could be a doublet or a triplet of 
events within a ~ 0.7 degree error circle within a day, 
for models i?3oo and Daoo- In model D^qq the spectral 
sub-cases D{1) (no pairs) and D{2) (full thermalization 
due to pairs) bracket the range of possibilities, the answer 
being probably closer to the latter. For z ~ 10 the signals 
would be larger, while for — lOOM© the signal is a 
factor ^3 — 4 times smaller, and even smaller for Mh — 
SOMq . For Mh = 3OM0 , a model B30 would be very hard 
or impossible to detect, especially against a background 
of assumed Pop. I/II GRBs; but even for this mass, 
the signal might be a doublet from a redshift z ~ 10 for 
_D(1)30 (note that Z?(2)3o is not appropriate since no pair 
formation is expected for this mass). 



B. Population III GRB rates 

For calculating the diffuse flux we need to know the 
burst rate as a function of redshift, which is very uncer- 
tain due to the lack of observations of confirmed Pop. Ill 
objects of any kind. Population III GRBs are both much 
rarer and located at higher redshifts than the usually 
considered GRBs from the second and subsequent gen- 
erations of stars (i.e. Pop. I/II GRBs [Ml 112]). If we 
assume that the Pop. Ill GRB rate traces the Pop. Ill 
star formation rate (SFR), the observed all-sky GRB rate 



can be parametrized as (1 -I- z)dNQ^^/dz = </'gkb(-^) = 
<l)'sFR(z)^GRBPph{z)dV/dz [1Q1IT4], where <^|%(z) is the 
Pop. Ill SFR per unit comoving volume (in units of 
Mq yr^^ Mpc~^), ecRB is the efficiency of the GRB 
formation (in units of Mq^), Pph{z) is the detection 
efficiency of photons for a specific instrument such as 
the Swift BAT, and dV/dz is the comoving volume ele- 
ment of the observed area per unit redshift. Adopting 
a specific result of the extended Press-Schechter simula- 
tion for 4>'sfr{^) ^-iid assuming that ecRB 10~^ 
(similar to that for the ordinary Pop. I/II cases) and 
Pph{z) ^ 0.3 for z ~ 10 — 20 (which corresponds to the 
case that the luminosity function of the Pop. Ill GRBs 
is similar to the ordinary Pop. I/II GRBs), one obtains 
(f>alfB{z) - 0.5 yr-i [H] (note that Swift BAT only cov- 
ers ~ 27r/3 of the sky, and then we use ecRB and ^q^^ (z) 
^ 6 times larger than those in [14] to obtain the all-sky 
GRB rate). 

However, for such bright bursts as we consider, the 
photon emission can be generally above the threshold 
of Swift BAT, i.e., Pphiz) - 1 [SmUj. Furthermore, the 
Pop. Ill stars can have a higher GRB formation efficiency 
eoRB- This may be written as ecRB = Vbeain^BHVenv, 
where rjbeam ~ is the jet beaming factor. In cases 
A and C, ribeam = 1/200 and in the cases B and D, 
Vbeam ~ 1/20000. The efficiency for the collapse to lead 
to a central black hole leading to a GRB is parametrized 
as esH (in units of M^^). Several simulations, e.g., [IH5] 
show that Pop. Ill stars are likely to have a top-heavy 
initial mass function (IMF, i.e. mass distribution), in- 
stead of a traditional negative index power-law (Salpeter) 
mass function. Stars in the mass range 140 ^ 26OM0 
undergo a disruptive pair-instability supernova explosion 
and leave no black hole remnant at all. Given the ma- 
jor uncertainties, we take here as a simple approxima- 
tive example a delta function IMF leading to black holes 
of Mh = 3OOM0, Mh = 100 A/q or 3OM0, so that the 
black hole formation efficiency is assumed to be as high 
as eBH ^ 0(1), 0(3) or o(10) per lO'^ solar masses. The 
remaining parameter rjenv denotes an efficiency factor re- 
lated to the environment under which the GRB jet can be 
formed. According to |14j , a requisite is that the envelope 
of the star be removed by a binary stellar companion in 
order to let the jet break out, which would suppress the 
GRB efficiency by a further order of magnitude. This 
condition implicitly assumes the jet durations of order 
30-100 s in the GRB frame known from low redshift ob- 
servations. In our model, however, the jets last td > lO^s, 
which is enough to break through even rather large stars, 
such as Pop. Ill without significant envelope mass loss. 
As an effective upper limit, we can assume r]env ~ o(l). 
The above represents an optimal theoretical case, which 
can be used as an upper limit, eoRB ^ 5 x 10~^ -^0^ 
for cases A300, C300, and ^ 5 x 10^^ ^^©^ for cases 
B300, Dsoo- Taking into account the larger Pph{z), we 
have a theoretical (and highly uncertain) upper limit of 
(f'GRB 10 yr~^ for cases B,D, and in principle two 
orders of magnitude higher for cases A, C. 
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We consider now some possible indirect observational 
constraints on Pop. Ill GRB rates, since direct obser- 
vational constraints are so far not available. The after- 
glow of the Pop. Ill GRB as well as its prompt emis- 
sion is bright enough to be detected by Swift and by 
ground based optical/near-IR telescopes |9l[T0]. On the 
other hand, the current observed GRB rate at z > 6 is 
~ 0.6 yr-i, i.e., 3 GRBs (GRB 050904,GRB 080913, and 
GRB 090423) during the 5-yr operation of Swift. This 
may be partly because the optical and near-IR observa- 
tions are more difficult at redshifts z > 6 due to the Lya 
absorption in the intergalactic medium. According to the 
current statistical data [53], only a small fraction ~ 25% 
of GRBs detected by BAT have redshift determinations, 
because of bad conditions for optical and near-IR obser- 
vations (not only the Lya drop-off effect but also condi- 
tions such as weather as well as dust extinction in the host 
and our galaxy) . Although the GRB redshift determina- 
tion rate may be a function of redshift (i.e., higher red- 
shift GRBs suffer stronger Lya drop-off), we may crudely 
estimate that the intrinsic GRB rate at z > 6 detected 
by BAT is - 0.6/0.25 = 2.4 yr'^. Since the BAT covers 
only ~ 27r/3 of the sky, the estimate of the GRB rate at 
z > 6 from the isotropic sky can be ~ 2.4 x 6 ~ 14 yr~^. 

Another indirect constraint on the rate of GRBs from 
the Pop. Ill very massive stars may be their production of 
very long and bright radio afterglows [lOj . The predicted 
durations at 1 GHz can be typically as long as 200 yr~^ 
for z > 10 radio afterglows arising from such GRBs. A 
comparison of the NVSS catalog (spanning over 1993- 
1996) and the FIRST catalog (spanning over 1994-2001), 
which effectively cover ~ 1/17 of the sky, did not find 
any radio transient sources which could be due to GRB 
afterglows with significant flux changes over timescales 
of ~ 5 yr [54l [55] . This indicates that the isotropic rate 
of Pop. Ill GRBs is constrained to be < 17/5 = 3.4 yr"^ 
From the above two rough estimates, we adopt in this pa- 
per a conservative observational constraint on the mas- 
sive Pop. Ill GRB rate of nf, < 3 yr^^. 



C. Diffuse neutrino flux and detectability 

With the above estimates of the rate of Pop. Ill GRBs 
we can now evaluate the diffuse neutrino flux expected 
over a given integration time, for comparison with the ca- 
pabilities of some current or planned neutrino telescopes 
such as IceCube [13], ANITA [SS] and ARIANNA [13]. 

In Figure 8 we plot the diffuse neutrino flux from 
Mh = 3OM0 Pop. Ill GRBs, the lowest mass case con- 
sidered, averaged over a period of a year of observation, 
in units of GeV cm~^ sr~^. This figure assumes, for 
illustrative purposes, a conservative GRB rate of = 1 
observed event per year (taking into account beaming ef- 
fects), and assumes that the bursts occur predominantly 
at a given redshift (or narrow range of redshifts) indi- 
cated in the figure, up to an upper limit of z = 70 for 
Pop. Ill formation [38] . The nominal case of one burst 



per year implies of course an anisotropic flux, but for 
a multi-year integration time, an averaged diffuse flux 
would be approximated. Also, if the rate were larger (e.g. 
in accord with the above observational indirect limit, say 
TT-b ~ 3yr~^), the diffuse fluxes would be higher by a fac- 
tor Tlf,. 

The neutrino spectra shown in Fig. 8 are for the 
cases i?3o and D(l)3o discussed above, since these are 
the preferred candidates for detection. However, for this 
Mh = 3OAf0 sub-case, even from z = 10 the predicted 
fluxes just approach the IceCube 5 year sensitivity [13] 
or the ARIANNA 5 year sensitivity [13] (if in the lat- 
ter we extrapolate the 6-month values to 5 years using a 
t^^^ scaling, which could be too conservative if the sen- 
sitivity is signal limited). For this = 3OM0 mass we 
have used the spectral case D{1) instead of D{2) because 
the lower luminosity leads to a 77 compactness too low 
to lead to significant thermalization. The diffuse fluxes 
for the cases A and C are not shown, since they are too 
low to be of observational interest. The relatively larger 
chance of detection in cases B and D is attributable to 
the smaller jet opening angle and/or a higher external 
medium density. 

In Figure 9 we plot the diffuse neutrino flux from 
Mfi — IOOMq Pop. HI GRBs, averaged over a period 
of a year of observation, assuming again a GRB rate of 
ni, = 1 event per year (including beaming effects) and 
assuming they arise from various redshifts. For this and 
the higher mass case we show the B and D{2) spectral 
models, the latter assuming that 77 effects have thermal- 
ized the photon spectrum. In this mass case, 5 years of 
observations would appear to make detection feasible if 
they arise predominantly from redshifts z ~ 10, or from 
z ~ 20 for a larger ~ 3 yr^^. 

Figure 10 shows the diffuse neutrino fluxes predicted 
for the highest mass case Mh — 3OOM0, averaged over 
a year of observation for a GRB rate of rib = 1 events 
per year (beaming effects included), from various red- 
shifts. The spectral cases shown are again the B and 
D{2) models, which give the highest fluxes (cases A and 
C giving observationally negligible fluxes, even at this 
high mass [0] ). For these high mass models, 5 years of 
observations would make detection feasible even if they 
arise predominantly from redshifts z ^ 20 and the rate 
is as low as n;, = 1 event per year. 

These Pop. Ill GRB diffuse neutrino fluxes of Figs. 8, 
9 and 10, for all three mass values and for ni, = 1 yi~^, do 
not exceed the Waxman-Bahcall cosmic ray limit (includ- 
ing evolution in redshift ^ , line marked "WB bound" ; 
see also [SS]) or the nominal GZK cosmogenic neutrino 
flux [iSWT] (line labeled "GZK"), and they are compati- 
ble with the observational constraints set by the ANITA- 



[0] In order to reach the minimum detectabihty with IceCube or 
ARIANNA, even for = SOOM© one would need for cases A 
and C a rate of at least 100 yr~^ if the redshift is as low as z = 7, 
and an even higher rate if z is higher. 
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2 mission [T21|5n]- The atmospheric neutrino background 
decreases extremely steep and is only relevant here at en- 
ergies < 10^-5 GeV, where its value is - 10"^ - 10"^ 
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Figure [8]: Diffuse neutrino flux in units of GeV 
cm~^ s^^sr^-'^ for GRBs of = 3OM0, averaged over 
a year in the observer frame, based on the nominal 
assumption of a rate n^, — \ yr^^, for the cases Sao find 
iD(l)30 (the latter because the compactness parameter 
is not large enough as to thermalize the photons to 
produce a -D(2)3o case). We have assumed different 
typical source redshifts. Also shown are the atmospheric 
neutrino background, the IceCube 5 year limits, the 
ARIANNA 6 month limits and the ANITA II 45 day 
limits. 
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Figure [9]: Diffuse neutrino flux in the same units 
as Fig. 8, for GRBs of — IGOM©, averaged over 
a year in the observer frame, based on the nominal 
assumption of a rate = 1 yr~^, averaged over a year 
in the observer frame, for the cases -Bioo and £'(2)ioo, at 
different typical source redshifts, with the atmospheric 
neutrino background, the IceCube 5 year limits, the 
ARIANNA 6 month limits and the ANITA II 45 day 
limits. 
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Figure[10]: Diffuse neutrino flux in units of GeV 
cm-2 s-^sr-i, for GRBs of My, = SOOM©, averaged 
over a year in the observer frame, based on the nominal 
assumption of a rate nf, = 1 yr~^, for the cases S300 and 
_D(2)3oo, at different typical source redshifts, with the 
atmospheric neutrino background, the IceCube 5 year 
limits, the ARIANNA 6 month limits and the ANITA II 
45 day limits. 



IV. DISCUSSION 

The earliest macroscopic objects to arise out the 
primeval Universe plasma arc the first generation Pop- 
ulation III stars. These are potentially invaluable probes 
of the cosmology in the epochs 2: ~ 10 — 30 leading to 
the currently detected Universe. The lack of heavy ele- 
ments at the formation of these first objects has led to 
expectation that they are very massive objects, in excess 
of hundreds of solar masses [1]-^ , although there is de- 
bate about the possibility of much lower masses [5] . The 
almost complete lack of observational data about these 
earliest objects is a major problem. However, objects 
more massive than about SOM© may collapse to become 
gamma-ray bursts (GRBs) based on what is observation- 
ally known at lower redshifts, and these could serve as 
invaluable tracers for this first generation of star forma- 
tion. These GRBs are expected to be bright 7-ray, X-ray 
and optical/infrared sources, but optical/IR detections 
providing redshift determinations are largely prevented 
by the Ly-a absorption from the intervening intergalac- 
tic gas, while gamma- or X-ray observations would be 
hard to disentangle from those of lower luminosity and 
lower redshift GRBs. For this reason, predictions about 
the expected high energy neutrino spectrum from Pop. 
Ill GRBs, providing a new channel completely free of 
absorption by the intervening medium, could provide an 
invaluable handle about the rate and characteristics of 
these objects, and through them, about the formation of 
the earliest stars and structures. 

One of the largest uncertainties is the initial stellar 
mass function. If the Pop. Ill star masses are as large 
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as often assumed [IHSl [S] the largest black hole masses 
Mh ~ IO^-^Mq used here are probably a conservative 
upper limit to the stellar core collapse remnants. The 
corresponding Pop. Ill GRB luminosities (eq. ([2])) used 
here are based on the electromagnetic extraction of the 
rotational energy of a fast rotating black hole, and the av- 
erage beaming angle 9j ^ 10~^ assumed is comparable to 
the values inferred for high luminosity Fermi bursts |60j . 
There is no agreement on the value of such parameters; 
for instance, ours differ from those of |61j who assumed a 
lower jet luminosity (i.e., af^/a/3 ~ 2 x 10~^) and larger 
jet angle dj ^ 10^^ calibrated on those of Pop. I/II low 
redshift GRBs, implying a larger energy input into the 
stellar envelope and less into the emergent jet. On the 
other hand, the higher jet luminosity and narrower jets 
adopted here lead to a higher ratio of emergent jet to 
envelope energy. 

Other uncertainties concern the rate of occurrence with 
redshift of Pop. Ill stars, the resulting Pop. Ill GRB 
rate, and the external medium density into which the 
GRB relativistic jet expands. For the formation rates 
we have based ourselves on recent numerical simulations 
and theoretical estimates (j |IIIB ), moderated by current 
upper limits on the number of possible unidentified high 
redshift bursts in the currently observed sample, as well 
as the possible contribution of their afterglows to existing 
transient radio source observations. A conservative up- 
per limit that we have considered is one Pop. Ill burst, 
on average, per year. The generally considered range of 
redshifts of occurrence of Pop. Ill stars is 10 < z < 30, 
with a likely typical redshift of z ^ 20. As far as the 
external gas densities, we have used I < n < lO^cm"'^, 
spanning a range of plausible values. Among the cases 
considered, those with high mass, high luminosity and 
narrow beam angles lead to high isotropic equivalent en- 
ergies, which together with moderate to high external 
densities result in relatively high fluxes. 

In view of the considerable uncertainty concerning the 
appropriate parameter values, the predicted fluxes must 
be considered as nominal values, to be tested via obser- 
vations in order to narrow down the range of possibilities. 

Our calculation of the proton acceleration and the pro- 
duction of high energy neutrinos via photo-meson in- 
teractions in the GRBs own photon field involve new 
features not present in previous calculations. Typical 
calculations of the neutrino emission from GRBs, e.g. 
[HIES mug, considered lower redshift Pop. I/II ob- 
jects (c.f. [64], and generally assumed acceleration in 
purely hydrodynamic internal shocks (c.f. |32|). This 
is also the kind of bursts assumed for the recent Ice- 
Gube observational upper limits |65) . A basic difference 
in our case is that Pop. Ill GRBs are thought to in- 
volve MHD jets leading to prominent photospheric and 
external shock emission [5Hl0j. where the target photon 
field is different. Also, due to the higher luminosity the 
photon density and pair formation effects are more im- 
portant. This, together with a more detailed treatment 
of the photo-meson cross section and multiplicity results 



in different neutrino spectra. 

The prospects for detecting neutrinos from single Pop. 
Ill GRBs with IceGube [II] (and KMSNeT [SS] or AR- 
lANNA [nil) appear realistic, provided they have large 
masses compared to their low redshift counterparts, and 
provided they are efficient proton accelerators. They are 
also limited the case where the external medium density 
encountered by the jets is high, n > 10^ cm^'^. Such 
values, although highly uncertain, are within the range 
of what is expected from numerical simulations. As dis- 
cussed in QUA up to once a year, at > PeV energies a 
3OQM0 GRB at z '-^ 20 could yield a doublet or a triplet 
of events over a time of a day in IceGube, and in some 
cases even a 3OM0 GRB at z ^ 10 could produce a dou- 
blet over a day. In the range 10-300 PeV these signals 
would be above a diffuse background from low-redshift 
Pop. I/II GRBs, and also of GZK cosmogenic neutri- 
nos of a different origin (while atmospheric neutrinos are 
not important at these energies). An accurate evalua- 
tion of the signal to noise ratio would however depend 
on model considerations for these backgrounds, which is 
beyond the purpose of this paper. While doublet and 
triplet searches have been done by IceGube over shorter 
timescales, < 100 s, recent extensions of such searches to 
multiplets over ~ day windows are of great interest for 
the signals discussed here. 

The detection of single sources would be aided if there 
were a simultaneous electromagnetic detection. The 
gamma-ray and X-ray flux would in principle be de- 
tectable [HI jTU], but for the long observer frame du- 
rations td,obs ^ 10^ s ^ 1 day the rise-time is very 
gradual and poses difficulties for normal triggering al- 
gorithms [To]. Optical detections are out of the ques- 
tion, due to the blocking by the intergalactic Ly-a ab- 
sorption [37' '68' . However, infrared L-band (3.4/im) or 
even K-band (2.2/im) detections may be possible for some 
models and redshifts. For instance, using the photon 
fluxes of Figs. 2 and 3, a rough estimate indicates that 
for Mh = 3OOM0 at z = 20, neither models i?3oo nor 
D(2)3oo are detectable in K but are detectable in L at 
the level of ~ 3 and 5, respectively, which is very 
bright. Scaling for an Mh — 30Mq burst, one would ex- 
pect at z = 20 neither B30 nor 0(1)30 to be detectable in 
K, but to be detectable in L at ^ 5 and 6 respectively. 
The L-band at z = 20 corresponds to source-frame UV 
frequencies, so there could be some intra-source absorp- 
tion making these dimmer. Observations of z > 8 GRBs 
[69] and galaxy candidates [70] do not show evidence for 
dust, although atomic or molecular resonant scattering 
could conceivably have some effect. Observations in the 
L-band, while not common, are being done with a few 
telescopes, but these generally have a small fleld of view, 
so one-day transients once a year such as described above 
could very easily have gone unnoticed. A detection at 
these wavebands would be dependent on having alerts 
based on automated triggers from gamma-ray. X-ray or 
neutrino signals. 

The prospects for detection of the diffuse neutrino flux 
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are also encouraging, under the above caveats of large 
masses, efficient proton acceleration and high external 
densities. Using a conservative Pop. Ill burst rate of 
nf, ~ 1 yr~^ and assuming black hole mass Mh ~ SOOMq 
or even M/,, ^ lOOAf© at z < 20 would lead to a diffuse 
flux, averaged over five years, which at energies > 1 PeV 
is within the reach of IceCube and ARIANNA, even in 
the presence of a diffuse neutrino background from lower- 
redshift Pop. I/II GRBs and GZK cosmogenic neutrinos 
of a different origin. The Pop. Ill GRB diffuse neutrino 
flux signals have a spectrum which differs significantly 
from that of the backgrounds mentioned above. Thus, 
if a sufficiently large number of events is accumulated, 
the spectrum should help to distinguish between these 
signals and the backgrounds. 

We note that the cosmic rays from these Pop. Ill GRB 
sources, after p7 losses within the sources and also in 
the CMB once outside of them, do not provide a signif- 
icant contribution to the diffuse cosmic ray background. 
Similarly, the neutrinos they produce do not contribute 
significantly to the GZK cosmogenic neutrino fraction. 

Thus, in five years or maybe less, IceCube would be 
able to rule out massive GRBs whose formation rate is 
rib 1 yr^i and Mu ^ SOOM© at redshifts z ~ 20, or 
nb ^ 3 yr""'^ and Mh ^ lOOAf© at z ~ 20, based on diffuse 
PcV neutrino measurements. The same measurements 
which are currently setting constraints on the previously 
considered Pop. I/II GRB diffuse background [12 [5T] 
will also be able to set constraints - or perhaps confirm 
- such models as considered here for the Pop. Ill GRBs 
and their environment. 

Since the neutrino detection of either individual Pop. 
Ill GRBs or their diffuse neutrino background is only 
possible for large black hole masses, implying progenitor 
stars of > 3Mh > IOOMq , large area neutrino ex- 
periments at energies > TeV would be able to address 
the currently unresolved question of whether Pop. Ill 



stars have very large masses or perhaps more modest 
masses approaching solar values. In the latter case, core 
collapse black holes and GRBs may not follow from the 
demise of this first generation of stars; the first black 
holes may arise from subsequent stellar generations, with 
smaller black hole masses, which would be much less lu- 
minous than those discussed here. (An ancillary impli- 
cation would be that a faster growth or coalescence rate 
would be needed to go from such later low mass black 
holes to the supermassive ones inferred in early quasars) . 
In the former case, the presence already in the Pop. Ill 
era of large mass ^ 30 — IOOMq black holes resulting in 
anomalously luminous GRB would provide a head-start 
for the growth into supermassive black holes, as well as 
provide information about the GRB physics in an ex- 
treme regime, testing questions of jet physics and prob- 
ing the star forming medium composition and density. 
Perhaps even more interestingly, PeV neutrino measure- 
ments could provide the first positive detections of the 
earliest massive structures to form in the Universe, at 
the so far unexplored 10 < z < 20 range of redshifts. 
Such measurements would be invaluable for a better un- 
derstanding of the earliest generation of stars, tracing the 
cosmic structure formation and the environment condi- 
tions at the dawn of the Universe. 
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